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ABSTRACT 

Newly born magnetars are promising sources for gravitational wave (GW) detec¬ 
tion due to their ultra-strong magnetic fields and high spin frequencies. Within the 
scenario of a growing tilt angle between the star’s spin and magnetic axis, due to the 
effect of internal viscosity, we obtain improved estimates of the stochastic gravitational 
wave backgrounds (SGWBs) from magnetic deformation of newly born magnetars. We 
find that the GW background spectra contributed by the magnetars with ultra-strong 
toroidal magnetic fields of 10^^ G could roughly be divided into four segments. Most 
notably, in contrast to the background spectra calculated by assuming constant tilt 
angles y = 7r/2, the background radiation above 1000 Hz are seriously suppressed. 
However, the background radiation at the frequency band ^ 100 — 1000 Hz are mod¬ 
erately enhanced, depending on the strengths of the dipole magnetic fields. We suggest 
that if all newly born magnetars indeed have toroidal magnetic fields of 10^^ G, the 
produced SGWBs should show sharp variations with the observed frequency at several 
tens to about 100 hertz. If these features could be observed through sophisticated de¬ 
tection of the SGWB using the proposed Einstein Telescope, it will provide us a direct 
evidence of the tilt angle evolutions and further some deep understandings about the 
properties of newly born magnetars. 
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1 INTRODUCTION 

Neutron stars (NSs) are one of the most promising tar¬ 
gets for gravitational wave (GW) detection with the ground 
based GW interferometers such as LIGO, Virgo, GEO600, 
advanced LIGO (aLIGO) and the proposed Einstein Tele¬ 
scope (ET; see Sathyaprakash & Schutz 2009 for a review). 
In general, GWs can be radiated from NSs because of many 
different reasons, e.g., magnetic deformation of the NSs 
(Bonazzola & Gourgonlhon 1996; Regimban & de Freitas 
Pacheco 2001, 2006a; Stella et al. 2005; Dall’Osso et al. 2009; 
Marassi et al. 2011), dynamical bar-mode instability (Lai 
& Shapiro 1995), r-mode instabilities due to GW radiation 
back-reaction (Andersson 1998; Friedman & Morsink 1998; 
Zhu et al. 2011), oscillations or flows excited by glitches 
(Abadie et al. 2011; Stopnitzky & Profumo 2014), collapse 
induced by phase transition to quark matter (Marranghello 
et al. 2002; Sigl 2006), and coalescences of compact binaries 
(i.e., NS-NS, NS-white dwarf, or NS-black hole; Schneider 
et al. 2001; Regimban & de Freitas Pacheco 2006b; Regim- 
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ban & Chauvineau 2007; Zhu et al. 2013). For a newly born 
NS of interest here, which if spins initially at the mass- 
shedding limit with a millisecond period and has an ultra¬ 
strong magnetic field, it could generate a prompt GW signal 
due to the bar- and r-mode instabilities and a continuous 
GW radiation due to a magnetically-induced quadrupole 
ellipticity. Such an ellipticity is mainly determined by the 
equation of state (EOS) of the NS, the magnetic configura¬ 
tion in the stellar interior and, most importantly, the mag¬ 
netic energy (Haskell et al. 2008; Dall’Osso et al. 2009). 

It is difficult, usually impossible, to detect GW radi¬ 
ation from magnetic deformation of the Galactic NSs. For 
canonical NSs such as the Grab pulsar and the central com¬ 
pact object harbouring in Gassiopeia A, their present mag¬ 
netic fields are obviously too low to sustain a sufficiently 
high ellipticity. Direct searches for GWs from nearby pul¬ 
sars has put an upper limit of e < 7 x 10~® at 95% confi¬ 
dence level for the pulsars (Abbott et al. 2010). Meanwhile, 
for the Galactic magnetars whose magnetic fields could be 
high enough, their long spin periods (P ~ 2 — 12 s; Olausen 
& Kaspi 2014) make them be below the optimal observa¬ 
tional bands of LIGO, Virgo and aLIGO. Nevertheless, for- 
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tunately, recent observations to superluminous supernovae 
(SLSNe) and gamma-ray bursts (GRBs) could open a com¬ 
pletely new window allowing us to search the GW signals 
from some extragalactic newly born NSs. The fittings to a 
remarkable number of SLSN light curves (Kasen & Bildsten 
2010; Inserra et al. 2013) and GRB X-ray afterglows owning 
a shallow decay/plateau phase (e.g., Dai & Lu 1998; Zhang 
& Meszaros 2001; Yu & Dai 2007; Yu et al. 2010; Dall’Osso 
et al. 2011; Rowlinson et al. 2013) robustly suggested that 
the NSs born there could be millisecond magentars (i.e., NSs 
with a high dipole magnetic field of ~ 10'^'^“'^® G). Such mil¬ 
lisecond magnetars can release a remarkable fraction of their 
rotational energy in a sufficiently short period to energize 
the SLSN and GRB outflows. 

On one hand, the millisecond periods of the newly born 
magnetars harbouring SLSNe and GRBs can determine the 
frequencies of their GW radiation to be within the interval 
of ~ 10^ — 10® Hz, which is just in the sensitive bands of 
LIGO and Virgo. On the other hand, in the sight of dy¬ 
namo models for magnetic field generation and amplifica¬ 
tion, an internal multipole (usually toroidal) magnetic field 
much higher than the surface dipole field can be simulta¬ 
neously formed due to the differential rotation of the newly 
born magnetars (Duncan & Thompson 1992; Cheng & Yu 
2014). The equipartition between the differential rotation 
and the toroidal field formation would determine the field 
strength to be on the order of ~ 10^® --gQ (Braithwaite 
2006). Such values can be supported by the observations to 
the giant flare event on 2004 December 27 from SGR 1806- 
20, which indicates an internal magnetic field of ~ 10^® G in 
that magnetar (Stella et al. 2005). The ultra-high internal 
toroidal magnetic field can lead to a very high quadrupole 
ellipticity, allowing the magnetars to produce strong GW 
radiation. In fact and specifically, the GW radiation of a 
newly born magnetar is sensitive to the tilt angle between 
its spin and magnetic axis and, as suggested by Dall’Osso et 
al. (2009), the tilt angle is actually time-dependent which is 
determined by the competition between the GW radiation 
and viscosity. Therefore, for an elaborate calculation of the 
GW radiation of newly born magnetars as well as their con¬ 
tributed stochastic gravitational wave background (SGWB), 
it is necessary to take the tilt angle evolution into account, 
in particular, with an ultra-high toroidal field. 

The SGWB due to magnetic deformation of newly 
born extragalactic magnetars has been widely investigated 
(Regimbau & de Freitas Pacheco 2001, 2006a; Marassi et al. 
2011), where however the tilt angle evolutions of magnetars 
are all neglected which could lead to some unrealistic conse¬ 
quences. Therefore, in this paper we revisit the calculation 
of the SGWB contributed by newly born magnetars by com¬ 
bining with a consideration of the tilt angle evolution. The 
paper is organized as follows. In Sections 2 and 3, the cal¬ 
culations of the GW spectrum from a single magnetar and 
the SGWB from all extragalactic newly born magnetars are 
introduced, respectively. In Section 4, we present our numer¬ 
ical results. Conclusion and discussions are given in Section 
5. 


^ An upper limit on the internal magnetic field can be set by the 
virial equipartition as = 1-4 x 

(Mastrano et al. 2011). 


2 GRAVITATIONAL WAVE RADIATION 

FROM MAGNETARS 

The ellipticity of magnetars has been investigated widely 
with different NS interior structures and different magnetic 
field configurations (Bonazzola & Gourgoulhon 1996; Cutler 
2002; Haskell et al. 2008; Ciolfi et al. 2009, 2010; Mastrano 
et al. 2011; Mastrano & Melatos 2012; Akgiin et al. 2013; 
Mastrano et al. 2013; Mastrano et al. 2015). Non-linear nu¬ 
merical simulations showed that the internal magnetic fields 
of a magnetar could probably have a ‘twisted-torus’ con¬ 
figuration (Braithwaite & Spruit 2004, 2006; Braithwaite 
2009), which is more complicated than the usually-assumed 
purely poloidal or toroidal structure. However, in any case, 
the toroidal field component is usually found to be the dom¬ 
inated one. Therefore, in this paper we simply take a pure 
toroidal magnetic field in the magnetar interior. In this case, 
the quadrupole ellipticity of the magnetar can be estimated 
as eb = —1.6 X 10“'^(i3t/10^® G)^, where Bt is volume- 
averaged strength of the toroidal field (Cutler 2002) and the 
minus sign indicates that the magnetar has a prolate shape. 

We follow the same procedure as Marassi et al. (2011) in 
deriving the GW energy spectrum emitted by a single newly 
born magnetar and the SGWB contributed by an ensemble 
of such magnetars. First, with a quadrupolar magnetic el¬ 
lipticity cb and tilt angle X: fhe rate of energy release from 
a newly born magnetar via GW radiation can be written as 
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^ ^ sm x(16sm x + cos x), 
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where I is the moment of inertia and ui the angular fre¬ 
quency. More specifically, the energy represented by the first 
and second term on the right side of the above equation 
would be mainly emitted into the GW frequency bands of 
a;/27r < i/e SS (jj/tt and i/g ^ u)I2'k, respectively. Therefore, 
the spectrum of the GW radiation can be approximated by 
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Due to the GW radiation and magnetic dipole radiation 
(MDR), the spin evolution of the magnetar can be deter¬ 
mined by 

2Ge%Iuj^ . 2 ■ 2 , 2 . 2 

^ = - ^ -sm x(16sm x + cos x) -sm x, 

where Bd is the strength of the surface dipole magnetic field 
at the magnetic pole and R the radius. 

Following Cutler & Jones (2001) and Dall’Osso et al. 
(2009), on one hand, the tilt angle can be decreased to zero 
due to GW radiation, no matter which shape the magnetar 
has. On the other hand, for a magnetar of prolate shape, 
the tilt angle is inclined to increase to tt/2 to minimize its 
precession energy. Therefore an increase of the tilt angle can 
be carried out through viscous damping of the precession 
energy with a time-scale (Dall’Osso et al. 2009) 


Td = 13.5 s 
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where Eb, P, and T represent the total magnetic energy 
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spin period, and temperature of the magnetar, respectively. 
In the following calculations, an analytical thermal history 
T{t) = 10® [t/1 yr + (10® is taken for modified 

Urea cooling of the magnetar (Owen et al. 1998). Then fol¬ 
lowing DairOsso et al. (2009), we have 


dsin X 
dt 


cos® X sin X 
sin® X 


20 * 2 j i ■ 

—smx. 


( 6 ) 


It should be noticed that the damping time-scale presented 
in Equation © is obtained with the bulk viscosity of neu¬ 
tron matter. In fact, as the temperature of the magnetar 
decreases to ~ 2 x 10® K (Page et al. 2004) at the time of 
~ 5 X 10® s, the neutron matter in the interior would enter 
into superfluid phase and thus the bulk viscosity disappears. 
Nevertheless, for temperatures lower than ~ 10® K (Charnel 
& Haensel 2008), a solid stellar crust can be formed to con¬ 
tribute a new damping effect due to the coupling between 
the core and crust (Alpar & Sauls 1988). The corresponding 
damping time-scale reads Tcc = nP/ee, where n represents 
the number of precession cycles (Jones 1976; Cutler 2002). 
It is difficult to determine the value of n in theory. For rela¬ 
tively slow rotations, n is estimated to be n 10®“^ (Alpar 
& Sauls 1988; Cutler 2002). In any case, for cb ~ 0.01, it 
could be reasonable to assume that tilt angle can be in¬ 
creased to 7r/2 immediately along with the crust formation. 

From the above equations, we can calculate the evolu¬ 
tions of the spin frequency i^ns(= a)/27r) and the tilt angle x 
of magnetars with different parameter values, as presented 
in Fig. [T] In principle, the spin-down of the magnetars is 
controlled by both MDR and GW radiation. To be specific, 
for Bt <C 10^^ G the MDR braking can always play a dom- 
inatlve role in the spin-down history, whereas for Bt = 10*^^ 
G the early spin evolution would be completely changed by 
the strong GW radiation. These results confirm the conclu¬ 
sions of DairOsso et al (2015). Of more interests, in agree¬ 
ment with DaU’Osso et al. (2009), our results show that the 
tilt angles of the magnetars for Bt <C 10^^ G can rapidly 
grow to 7r/2 only in a few seconds, because of the short vis¬ 
cous damping time-scale for a millisecond spin period. For 
Bt = 10*^^ G, even though the viscous damping is suppressed 
by the strong GW radiation, the tilt angle can still be in¬ 
creased to 7r/2 at about 5 x 10® s due to the curst formation 
and the consequent core-crust coupling. In other words, for 
newly born millisecond magnetars, the orthogonal configu¬ 
rations with X = 7r/2 can be built quickly In all situations, 
which are the most benefit for their GW radiation. 

With the above stellar evolutions, we show the corre¬ 
sponding GW radiation spectra in Fig. [21 where some re¬ 
sults for constant tilt angles x = ^/2 are also presented for 
comparisons. As shown, for Bt <C lO^’^ G, the influence of 
the tilt angle evolution on the GW radiation spectrum is 
very limited, since the angle changes too quickly. However, 
for Bt = G, some significant changes in the spectrum 
appear, in particular, in the high-frequency domain. To be 
specific, according to the contributor of radiation, the GW 
spectrum can be roughly divided into four segments. For in¬ 
stance, in the case of Bd = G and Bt = 10^^ G, follow¬ 
ing Equations ([2]) and Q, and defining A = 



Figure 1. Evolutions of the spin frequency i^NS (upper panel) and 
tilt angle x (lower panel) of newly born magnetars. The strengths 
of the surface dipole and the internal toroidal magnetic fields are 
shown in the legend. The initial conditions for the magnetars are 
Pi = 1 ms, Xi = l°i and Ti = 10^® K. Other quantities for the 
magnetars are taken as M = 1.4 Mq, R = 10® cm, and I = 10^® 
g cm^. 


C = we have 

Asin®x[R -b <7(1 -b 15sin®x)]“\ • • • 

1000 < r-e < 2000 Hz 
Acos®x[R -b 4(7(1 + 15sin®x)]~^, • • • 

222 < t-e < 1000 Hz (7) 

Acos®x[R -b 4C'(1 + 15sin®x)]”^ ■ • ■ 

+ A(B + 16(7)"\ 111 <Ue< 222 Hz 
A{B + 16Cy^, r/e < 111 Hz. 

From the above equation we know that the first and fourth 
segments of the spectrum are contributed by the radiation 
at 2z/ns- While the second segment Is dominated by the ra¬ 
diation at i/NS and the third part contains the contributions 
from iZNS and 2:/ns- On the other hand, the GW energy 
spectrum calculated for x = 7r/2 can be expressed as 

= A(B + 16(7)"\ iZe < 2000 Hz. 
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Figure 2. The GW energy spectra emitted by newly born mag- 
netars with different strengths of the dipole and toroidal magnetic 
fields. The red and green lines represent the results calculated by 
involving the tilt angle evolutions. The dipole and toroidal fields 
are the same as those in Fig. ^ The grey and blue lines show the 
spectra derived only by assuming constant angles x — but 

with the same strengths of magnetic fields as in the cases of x 
evolution, as shown in the legend. 


Dividing the four segments of Equation 0 with Q, 
respectively, one can see the GW radiation is reduced by 
sin^X — Ibsin^x (due to B <C C) at the frequency 
band 1000 < Ve < 2000 Hz; enhanced by cos^x — 4 

(due to B <C C) at 222 < r'e < 1000 Hz; enhanced by 
4-4 < 1 + -l^cos^x < 5 (due to B < G ^ B ~ C) at 
111 < iZe < 222 Hz; unchanged at Vg < 111 Hz as compared 
with the result derived by assuming x = 7r/2. These analytic 
results accurately account for the enhancement in the GW 
energy spectrum as shown in Fig. [21 

Likewise, in the case of Ba = 10^® G and Bt = 10^^ 
G, the GW energy spectrum could also be divided into four 
segments, namely, 1000 < Vg < 2000 Hz, 140 < iZe < 1000 
Hz, 70 < iZe < 140 Hz and Vg < 70 Hz. From high- to low- 
frequency band the expression for each segment is the same 
as Equation 0. Owing to the rather small tilt angle at early 
period, the GW radiation is also suppressed at 1000 < Vg < 
2000 Hz. However, as Ba increases, much more rotational 
energy is released into MDR. Hence, the enhancement in 
the GW radiation is not as remarkable as in the previous 
case. Quantitatively, the GW energy spectrum is enhanced 
by 1 < ^^cos^x < 3.4 (B < G B > G) and 1 -b 
cos^x ~ 2 (B 3> G) at 140 < Vg < 1000 Hz and 
70 < I'e 140 Hz, respectively. Moreover, we confirmed the 
result of Marassi et al. (2011), which suggested that at low 
frequency band where MDR dominates the spin-down, the 
GW energy spectrum is reduced by two order of magnitudes 
while increasing Ba from 10^^ to 10^® G. 


3 STOCHASTIC GW BACKGROUND 

Generally, the SGWB is denoted by a dimensionless quan¬ 
tity, flGw(r'obs), which actually represents how the GW en¬ 
ergy density is distributed with the frequency t'obs in the 
observer frame. Following Ferrari et al. (1999), this quantity 


takes the form 

DGw(r'obs) ~ ^obsFi^ (i^obs), (9) 

(rpg 

where pc = SHq/StvG is the critical energy density needed to 
close the universe. The GW flux at the observed frequency 

^obs is 

FvCohs) = j f^{i'ohs)dR{z), (10) 

where fvCohs) is the GW energy flux per unit frequency of a 
single source, dR{z) is the magnetar formation rate located 
between z and z + dz. The GW energy flux emitted by an 
single magnetar is related to its GW energy spectrum as 
(Regimbau & de Freitas Pacheco 2006a) 

where di is the luminosity distance, and the GW frequency 
in the source frame has the form i/g = r'oba(l + z). 

The number of NSs formed per unit time out to redshift 
2 can be estimated by integrating the cosmic star forma¬ 
tion rate (CSFR) density, p*(z), over the comoving volume 
element, and taking into account the restriction associated 
with the initial mass function (IMF). Of all the NSs, 8—10% 
are considered to be magnetars (Bd > 10^"* G) according to 
the results of Regimbau & de Freitas Pacheco (2001) and 
Popov et al. (2009), who derived this ratio through popula¬ 
tion synthesis methods. In this paper, we assume the ratio 
to be ^ = 10%. As a consequence, the number of magnetars 
formed per unit time out to redshift z within the comoving 
volume can be written as 

fZ i-w rznmax 

RO = i / P*C) — dz / $(m)dm, (12) 

h dz' 


where dV/dz is the comoving volume element, and 4>(m) 
is the IMF. We take the CSFR density model suggested in 
Hopkins & Beacom (2006), which has refined the previous 
models up to redshift z ~ 6, and derived a parametric fit 
expression for the CSFR based on the new measurements of 
the galaxy luminosity function in the UV and far-infrared 
wavelengths. The CSFR density can be expressed as 


p*(z) = h 


0.017-bO.lSz 
l + (z/3.3)5-3 


Moyr ^Mpc 


(13) 


where h = 0.7. 

The comoving volume element in Equation m takes 
the following form (Regimbau & de Freitas Pacheco 2001, 
2006a; Zhu et al. 2011) 


dV , c r(z) 

- — 471 -- ^ - 

dz Ho E{n, z) ’ 


where B(D, z) = ■\/Dm(l + z)3 -|- Da, and r(z) = ci_L/(l -b z) 
is the comoving distance. Ho is the Hubble constant. We take 
the ACDM cosmological model with Ho = 70 km s“^Mpc“^, 
Dm = 0.3 and Da = 0.7. The standard Salpeter IMF 
4>(m) = is adopted with x = 1.35, where A 

is a normalization constant, determined by the relation 
m$(m)cim = 1. It should be noted that the mass 
range of magnetar progenitor is still controversial (see Fer- 
rario & Wickramasinghe 2008; Davies et al. 2009). However, 
we mainly focus on the effect of the tilt angle evolution on 
the resultant SGWB. For simplicity, we take the same mass 
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range mmin = 8M©, rrimax = 40M© as Marassi et al. (2011) 
for the magnetar progenitors. 

Combining the Equations ®, Cni, OH), OH and d, 
one can obtain the SGWB contributed by the magnetically- 
induced deformation of an ensemble of newly born magne¬ 
tars 


flGw(r'obs) 


STvGvghs /"""PP p^,{z) dEg^ 

Jo {1 + z)E{Q.,z) due 


X 


pn 

J m. 




(15) 


The upper limit of the redshift integration is determined by 
the maximal GW frequency in the source frame and the 
maximal redshift 2 , of the CSFR model. To be specific, 
2:upp =min(2*, r'e.max/i^obs — !)• In deriving Equation (fT^ . 
the time-dilation effect has been involved, so the term (l-fz) 
appears in the denominator. 


4 RESULTS 

In Fig. 0] we plot the dimensionless energy density flew 
versus the observed frequency i/oba for newly born magne¬ 
tars with different strengths of toroidal and dipole magnetic 
fields. Since for a magnetar with Bt 10^^ G, its tilt angle 
evolution could not change the emitted GW energy spec¬ 
trum, the background spectrum contributed by an ensemble 
of such magnetars is also unchanged even their angle evo¬ 
lution is involved. However, for newly born magnetars with 
Bt — 10^^ G, the background spectra are obviously changed 
by the angle evolutions. These magnetars could keep rather 
small tilt angles for ~ 5 x 10® s, during which they have 
largely spun down. As a consequence, the GW radiation 
from 2 i/ns are seriously suppressed, resulting in the rather 
weak backgrounds above 1000 Hz. The background spec¬ 
tra show two cutoffs with frequencies at 1000 and 2000 Hz, 
respectively. Similar cutoffs are also presented in the spec¬ 
trum contributed by magnetars with constant tilt angles 
7r/60 (Marassi et al. 2011). Moreover, the background spec¬ 
tra contributed by newly born magnetars with Bt = 10^^ G 
can also be divided into four segments. 

Specifically, for Bd = 10^^ G and Bt = 10^^ G, the 
first segment of the spectrum is at the frequency band 
1000 < Vobs < 2000 Hz, the dimensionless energy density 
Hgw is reduced by a factor IGsin^y compared with result 
of X = tt/2. The evolution of Hgw with Vahs above 1000 Hz 
is mainly due to the tilt angle evolution at early periods. 
Initially, the tiny angles lead to rather small Hgw around 
2000 Hz. Subsequently, as the tilt angles increase and then 
keep almost constant values of 3°, Hgw increases and forms 
a plateau at 1000 < iZobs ^ 1250 Hz. However, such an evo¬ 
lution behaviour of the background spectrum is unlikely to 
be detected currently. The second and third segments of the 
background spectrum are at the band 222 < Vobs < 1000 Hz 
and 111 < Vobs 222 Hz, respectively. Quite the contrary, 
by involving the tilt angle evolution, Hgw at the two bands 
are respectively increased by about 4, and 4.4 — 5 times, 
following the analysis in Section 2. The enhanced SGWB 
could be detected by aLIGO at a few hecto-hertz. The main 
reason that leads to the enhancement is though the tiny tilt 
angle can suppress the GW radiation at 2iZNs of a single 
magnetar, more rotational energy is released into GW with 



Figure 3. Dimensionless GW energy density Ogw versus the 
observed frequency Vobs for newly born magnetars with different 
strengths of the dipole and toroidal magnetic fields, which are the 
same as in Figs.[T]and[2](see the legends in these figures). The red 
and green lines represent the SGWBs calculated by involving the 
tilt angle evolutions. The grey and blue lines show the SGWBs 
derived only by assuming x = rr/2, but with the same strengths of 
magnetic fields as in the cases of x evolution. Please note that the 
background spectra contributed by magnetars with Bd = 10^^ G 
and Bt = G are not shown for both the y evolution case 

and the y = 7r/2 case because they are too weak. The detection 
thresholds of LIGO (black solid line), VRIGO (black dashed line), 
aLIGO (black dotted line), and the proposed ET (black dash- 
dotted line), calculated using Equation (136) in Sathyaprakash & 
Schutz (2009) and assuming 1 yr observation time are also shown. 


emitted frequency at Ubis instead. The fourth segment is at 
iZobs 111 Hz, which has no differences with that derived 
by assuming y = n/2. 

The background spectrum produced by magnetars with 
Bd = 10*^® G and Bt = 10^^ G also contains four seg¬ 
ments, which are respectively 1000 < iZobs < 2000 Hz, 
140 < zzobs < 1000 Hz, 70 < zzobs 140 Hz, and Uobs 
Hz. Gompared with the result of y = tr/2, the first seg¬ 
ment is also reduced by IGsin^y, while the fourth segment 
keeps the same. Nevertheless, the second and third parts 
are enhanced by < 3 times and 2 times, respectively. Over¬ 
all, stronger Bd could lead to much more rotational energy 
release in the form of MDR, rather than GW radiation, fur¬ 
ther result in a relatively weak SGWB. As shown in Fig. 
O the background produced by newly born magnetars with 
Bd = 10^® G and Bt = 10^^ G may only be detected by 
ET. Actually, if the GW radiation braking is dominative 
in the early spin evolution of the magnetars, the resultant 
background spectrum is relatively flat at a few tens to a few 
hundred hertz, because dEg^/dUg oc Ve (GW radiation dom¬ 
inated) as compared to dEgw/dve oc zz® (MDR dominated) 
(Marassi et al. 2011). 

One can expect that if all newly born magnetars in¬ 
deed have ultra-strong toroidal magnetic fields of 10*^^ G, 
the resultant dimensionless GW energy densities would show 
sharp variations with the observed frequency at several tens 
to about 100 hertz. We suggest that if these features could 
be observed through precise detection of the SGWB by ET, 
it will provide us a key evidence about the early evolu¬ 
tion of tilt angles. In addition, this would provide us some 
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significant information about the properties of newly born 
magnetars, snch as their braking, cooling mechanisms and 
the viscosity of nnclear matter. Of conrse, the current non¬ 
detection of the SGWB by aLIGO mainly has the following 
reasons: (i) not all newly born magnetars have nltra-strong 
toroidal magnetic fields of 10^^ G; (ii) of all the newly born 
NSs, magnetars may have a ratio less than 10%; (iii) if all 
newly born magnetars have Bt — 10^^ G, then most of them 
should have Ba = G at least after their births. 


5 CONCLUSION AND DISCUSSIONS 

We have revisited the SGWB contributed by the 
magnetically-induced deformation of newly born magnetars 
involving the tilt angle evolutions based on the work of 
DairOsso et al. (2009). We find that for the magnetar with 
toroidal magnetic fields Bt <C 10^^ G, its tilt angle evolu¬ 
tion have no influence on the emitted GW energy spectrum. 
Hence, the SGWB from an ensemble of such magnetars is 
also unaffected by the angle evolution. However, for magne¬ 
tars with Bt = 10^^ G, their tilt angle evolutions could ob¬ 
viously change the GW energy spectra and the background 
spectra, which both could roughly be divided into four seg¬ 
ments. Specifically, the background radiation above 1000 Hz 
are suppressed by a factor of lOsin^y. The intermediate two 
segments range from 100— 1000 Hz, the background radi¬ 
ation of which are enhanced by about four times, and 4.4 — 5 
times for Bd = 10^^ G (< 3 times, and 2 times for Bd = 10^® 
G), respectively. Since more rotational energy of the mag¬ 
netars are released into MDR for higher Bd, leading to the 
relatively weak enhancement in the GW background radia¬ 
tion. The last segments are at iZobs < 100 Hz, the spectra 
of which are the same as the results of x = tt/2. As a con¬ 
sequence, involving the tilt angle evolutions could result in 
some more stronger and detectable SGWBs, if the newly 
born magnetars indeed have ultra-strong toroidal magnetic 
fields. We expect that the evidence of tilt angle evolutions of 
newly born magnetars may possibly be found in their SG¬ 
WBs. Moreover, sophisticated detection of the SGWB at the 
band below a few hundred hertz may put some constrains 
on the properties of these magnetars. 

Finally, two aspects that may affect the tilt angle evo¬ 
lution should be considered in detail in future work. As we 
know, the bulk viscosity of quark matter is much larger than 
that of neutrons, protons and electrons nuclear matter. As a 
result, if the magnetars are quark stars or even hybrid stars, 
their tilt angle evolutions should be different from those of 
classical NSs considered in this paper. The resultant SGWBs 
are probably also different. On the other hand, the elliptic- 
ity and magnetic energy of the magnetar with ‘twisted-torus’ 
magnetic field configuration are different from that derived 
based on the toroidal-dominated configuration (Cutler 2002; 
Haskell et al. 2008; Ciolfi et al. 2009, 2010; Mastrano et al. 
2011; Mastrano & Melatos 2012; Akgiin et al. 2013; Mas¬ 
trano et al. 2013; Mastrano et al. 2015). Hence, if the inter¬ 
nal fields of magnetars indeed have a ‘twisted-torus’ shape, 
we may expect different tilt angle evolutions, which further 
result in various SGWBs. 
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